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1 Introduction 



Linear Control systems are typically given by a differential equation, linking the state x and 
the control u 

x{t) = Ax{t) + Bu{t), t £ M>o, 

usually completed by an initial condition x(0+) = xq, and an algebraic equation linking state, 
control and the observation y 

y{t) = Cx{t) + Du{t) t G M>o, 

where A, B, C and D are matrices of appropriate sizes. Rewriting these equations as one 
system acting on the whole real line R instead of the positive half-line M>o we end up with 
an differential-algebraic system of the form 




where Oq denotes the derivative with respect to time and the initial condition for the state 
variable x transforms into an additional Dirac-(5-source term on the right hand side. Sys- 
tems of this form have been studied in the finite- and infinite-dimensional case in various 
works. In the infinite-dimensional case the operators B, C and D, acting on some suitable 
Banach- or Hilbert-spaces, are usually assumed to be bounded, while the operator A is a 
generator of a Co-semigroup. In this case the solution theory is rather straightforward. How- 
ever, in case of boundary control and observation it turns out that the operators B and C 
are in general unbounded and hence, a more sophisticated theory is needed. The classical 
approach is to consider so-called admissible operators B and C as it was done for instance in 
[ISKIllElElElEllITKIHllinill]. We focus on a class of hnear control problems, where the 
operators B and C are unbounded, but in our approach the admissibility for these operators 
has not to be verified. 

The solution theory provided in this article is based on the theory of evolutionary equations 
as they were considered in [9] . As it was shown in |12[ [TT] , linear control systems (including 
the case of unbounded operators B and C) are just a subclass of evolutionary equations. In 
this note we will generalize the solution theory presented in |12j to a broader class of so-called 
(linear) material laws. This generalization allows us to study control problems including delay 
terms. As in |12) we introduce abstract boundary data spaces, which enable us to formulate 
boundary control and observation equations without strong smoothness assumptions on the 
boundary of the underlying domain. Indeed, it will suffice to guarantee a Poincare-type esti- 
mate for the involved differential operators. 

Section [2] recalls some preliminaries on evolutionary equations, linear material laws and ex- 
trapolation spaces (so-called Sobolev chains) and we refer to [9l [TOl [6] for the proofs and a 
deeper study of the related topics. 

In Section [3] we introduce the notion of linear control systems, which will be a special case of 
the broader class of abstract evolutionary equations. In contrast to |12| we will generalize the 
class of possible control problems to the case of arbitrary material laws (while in |12| just the 
so-called (P)-degenerate case, cf. [^, was treated). We provide a well-posedness result for this 
class, which is in essence just an application of |9l Solution theory], and show the causality of 
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2 Preliminaries 



the solution operator. 

Boundary control problems are introduced in Section H] and we show that they fit into the 
abstract class of linear control problems introduced previously. In order to formulate bound- 
ary control and observation equations, without imposing strong smoothness constraints on 
the domain, we introduce abstract traces and recall the notion of abstract boundary data 
spaces. Finally, we apply our findings to a boundary control problem for the equations of 
visco-elasticity. 

2 Preliminaries 

In this section we recall the notion of evolutionary equations. Following [6] we begin to intro- 
duce the time derivative do as a boundedly invertible operator on an exponentially weighted 
L2-space. 

Definition 2.1. For ly G R>o we denote by Hi,^q{M) the space of all square-integrable func- 
tion^ with respect to the exponentially weighted Lebesgue-measure exp(— 2z^t) dt, endowed 
with the inner product given by 

im.fi ■■= J f{trg{t)eM-^'^t) dt if, 9 G ^.,oW). 
We define the operator do^u on Hiy^Q(R) as the closure of the derivative 

This operator is normal with 9^e5o,i/ = and hence, G Qido^^) with H^f^'i'/ll ^ ^ • If the 
choice of i/ is clear from the context we will omit the additional index i'. 

Clearly, we can extend the operator do to the space of //-valued functions HyfiiM; H), where 
is a Hilbert space, in a canonical way. 

Remark 2.2. 

(a) For u G H,yfi{M.) the function d^^u is given by 

X 

{dQ^u) (x) = j u{t) dt (j; e R a.e.). 

— oo 

This especially yields the causality of 9(7^0 

^Throughout we identify the equivalence classes induced by the equality almost everywhere with their repre- 
sentatives. 

■^A mapping F : //y,o(R; H) — > Hi^fi(9.\ H), where H is an arbitrary Hilbert space, is called causal, if for each 
a G R it holds X(-aa,a]{m)Fx(-oo,a]{'m) = X(-aa,a]{m)F, where by X(-oo,a]("^) we denote the operator on 
H^^o(M.; H) mapping a function / to the truncated function t i— >■ X(-oo, a] 
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(b) Let T : L2(M) — >■ ^2(1^) denote the Fourier-transform and define for ly > the operator 
e"'"™ : H^fl{R) L2(M) by (e"'""/) {x) = e"''^/(x) for almost every x e R, which is ob- 
viously unitary. Then we define the Fourier- Laplace-transform := Te~^^ : H^^q{]R.) — t- 
L2(1R), which gives tlie following spectral representation of the derivative do^i/'. 

do,u = >C*(im + u), 

where m denotes the "multiplication- by-the-argument" operator on -L2(M) with maximal 
domain. 

With the spectral representation of do we can define so-called linear material laws (cf. [9]). 

Definition 2.3. Let r > 0, H an arbitrary Hilbert space and M : Bc{r,r) — )• L{H) be 
bounded and analytic. Then we define the operator M ^ im+t/ ^ L2{R) for 1/ > ^ by 

and the linear material law M{dQ^) € L(H,yfi{R; H)) for > ^ by 

' \im -|- z^y 

Remark 2.4. Due to the analyticity of M we obtain by a Paley- Wiener result that the operator 
M{dQ^) is causal. 

Note that any densely defined closed linear operator A defined in a Hilbert space H gives rise 
to densely defined closed linear operator in Hi^fi{R; H) defined as the canonical extension of 
the operator acting as {Af){t) := Af{t) for all t G M and simple functions / taking values in 
the domain of A. Henceforth, we will identify A with its extension without further notice. 

Theorem 2.5 (Solution theory for evolutionary equations [9l Solution theory]). Let H be 

a Hilbert space and A : D{A) <^ H ^ H be skew-self adjoint. Furthermore let r > and 
M : Bc{r,r) — t- L{H) be analytic, bounded and assume that there exists c > such that for 
all z G B£{r, r) 

mez~^M{z)>c. (1) 
Then there exists z/q > such that for all u '>vq the evolutionary equation 



doM{d^') + Aju = f, (2) 

admits for every f £ Hj^^q(R\ H) a unique solution u G Hyfi{R\ H), which depends continuously 

on f. More precisely, G ^ ^doM^dQ^) -\- A^^ and the solution operator (^doM{dQ^) + A^ 
is causal. 

Next we introduce the concept of Sobolev-chains. For the proofs and further details we refer 
to [ini Chapter 2]. 

Definition 2.6. Let H he a. Hilbert space and C : D[C) <^ H ^ H he a densely defined, 
closed linear operator with G q{C). For k & 1^ we define Hk{C) as the completion of the 
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3 Abstract linear Control Systems 

domain D{C^) with respect to the norm \C^ ■ \h. Then {Hk{C))f^^^ is a sequence of Hilbert 
spaces with H].{C) i7fc_i(C) for k £ 7j. The sequence {Hi.(C))kez is called the Sobolev- 
chain of C . For each A; G Z the operator C : i?|^.|_(_i(C) C Hk+i{C) — )• Hk{C) possesses a 
unitary extension to Hf;^i{C), which will be again denoted by C. Furthermore H]f_{C)* can 
be identified with H-k{C*) for each k £ Z via a unitary operator. 

Remark 2.7. 

(a) Let Hq, Hi be two Hilbert spaces over the same field and A : D{A) C Hq — )■ Hi be densely 
defined, closed and linear. For each /c G Z the operator 

A : i^|fcn_i(|A| + i) C + i) ^ Hk{\A*\ + i) 

has a unique continuous extension to + i)- 

(b) Let (HkiC)) 

fcgz ^ Sobolev-chain associated to some operator C and A : Hi{C) — )■ H 
be linear and bounded, where H denotes an arbitrary Hilbert-space. Then the operator 
A'^ : H ^ H-i{C*) is defined as the dual operator of A, where we identify the dual of H 
with H and Hi{C)* is identified with i?_i(C*). 



3 Abstract linear Control Systems 

In this section we introduce the shape of linear control systems and show that they fit into 
the class of evolutionary equations introduced in the previous section. We consider a densely 
defined closed linear operator F : D{F) C Hq — t- Hi for two Hilbert spaces Hq and Hi. 
Furthermore let U and Y be Hilbert spaces, which will serve as control and observation space, 
respectively. 

Definition 3.1. Let Mi,i2 G L{Y;Hi) and Mi,2i G L{Hi;Y) for i G {0,1} and Mi,22 G 
L{Y;Y). Let ' ' 








Mi,02 ' 



M{z):=\ \0J \+z\ \Mi,i2j I {z€Bc{r,r)), (3) 

A/l,20 A'/i,21 ) Mi,22 

where K : Bc(r,r) — )• L{Hq © Hi) is a linear material law. An abstract linear control system 
Cm,f,b is an evolutionary equation of the form 

F*0\\ /x\ 
doM{dQ^)+\-F oil Uj=f + Bu. 

Here / G -ffjy^_oo(IR; Hq © Hi © Y) is an arbitrary source term, u G Hiy^Q{M.; U) is the control 
and B G L{Hy^Q{R; U); H,^^o{R; {Hq © i^i © Y))) is the control operator. We call an abstract 
linear control system well-posed, if the operator 

F* 0\ 

dQM{do^)+ I -F QH^,q{R;Ho®Hi®Y)(BH^,o(.^;Hq®Hi®Y) 
/ 
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possesses a densely defined bounded inverse for sufficiently large u. The continuation to the 
space H^fi{E.\ Hq © Hi © Y) of the inverse is called solution operator. 



It is clear that an abstract linear control system Cm,f,b is of the form ^ given in Theorem 
/ F* 0\ 

with A = { —F I . Hence, the solution theory for evolutionary equations is applicable. 
\ 0/ 

It is obvious that if M satisfies the condition ([1]), then so does N and the operator lHeMi^22 
is strictly positive definite. However, the latter is not a sufficient condition for the positive 
definiteness of 9\e z~^M(z). We define the operator J G L(Y; Hq © Hi) by 

^ 1 /Mi,02 + Mf20 

2 VA^i,i2 + M*2i 

Theorem 3.2. Let Cm,f,b be an abstract linear control system. Assume that D\ez-^N{z) > 
Co > and 9'leMi^22 > ci > 0. Assume that there is 6 > such that cq — S\\J\\ > and 
ci — ^\\J\\ > 0. Then C^[^f,b is well-posed and the solution operator is causal. 



Proof. Since, 



9\c 







Mi,02 

Mi,20 Mi,21 ) Mi,22 



we get for w := {x, (,,y) e Hq ® Hi ®Y 
D\z{z-^M{z)w\w) =yie(^z~^N{z) 

> Co 

>{co-m\\) 

The assertion then follows by Theorem 12.5 



J 
J* IHe Ml 22 



+ (9^eMi,222/|y) + 21He( Jy 



+ ci\y\'-2 

2 



+ Ci 



\J\ 



□ 



4 Boundary Control Systems 

This section is devoted to the study of boundary control systems. At first we show how 
boundary control and observation equations can be handled within the framework presented 
in the previous sections. This will mainly be done by a particular choice for the unbounded 
operator F. As it was pointed out in |12| Subsection 5.1], the resulting class of control systems 
can be interpreted as a generalization of a subclass of port-Hamiltonian systems (cf. [5]) to 
the higher dimensional case. Moreover, we recall the notion of so-called boundary data spaces, 
introduced in |12| Subsection 5.2], as well as abstract traces, which enable us to treat boundary 
values as suitable distributions belonging to some extrapolation space. 

First let us fix some notation. For Hilbert spaces Hq, . . . , Hn we define for i E {0, . . . , n} the 
operator 

TTHr-HQ®...®Hn^Hi 
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4 Boundary Control Systems 



as the orthogonal projection on Hi. Note that then 7r|^, is the canonical embedding from Hi 
to Ho® ...®Hn. 



We begin this section with an illustrative example. 

Example 4.1. Let C be an arbitrary domain and define the operators grad and div as 
the closures of 

gradlc-(n) : ^ L2{n) ^ ^2(0)" 

^ (5i</')i6{l,...,n} 

and 

divic^oo(n)" : c^inr c L2{nr ^ L2{n) 

n 
i=l 

respectively. These operators are formally skew-adjoint, i.e., grad C — ^div^ =: grad and 

div C — ^grad^ =: div . Then, using the extrapolation spaces of the operators | div | + i and 
I grad I + i we define the Dirichlet-trace and the Neumann-trace by 

7grad : ^^i(| grad | + i) ^ i^-id div | + i) 
u I— 7- ^grad — grad^ u 

and 

7div : ^^i(| div I + i) ^ F_i(| grad | + i) 
C I— 5- (div— div) 



respectively. Note that in the case of a smooth boundary, the distributions 7grad'U and 7divC 
for u G Hi{\ grad | + i) and C € Hi{\ div | + i) are supported on dil. More precisely with the 
help of the divergence theorem. 



(7gradw|C> = J u*C-n dS = (u|7divC>, 
dn 

where n denotes the unit outward normal and S the surface measure on (90. Note that 
Tgrad^^ = if and only if ti G D(grad) and 7divC = if and only if C G -D(div). 



In the rest of this subsection we generalize the concepts illustrated in the example above. For 
that purpose let Hq and Hi be two complex Hilbert spaces, G C Hq © Hi and D C Hi © Hq 
two densely defined closed linear operators, which are formally skew-adjoint. We define G ■= 
- (dY and D — -(g) 
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Definition 4.2 (Abstract traces). We define the abstract traces jg a-^d Jd by 

7G:ii'i(|G|+i)^F_i(|D|+i) 
(g-G^v 

and 

7D:^i(|i^|+i)-^i^-i(|G|+i) 
w H> ( D - d] w. 



Furthermore we define the abstract trace spaces as the image spaces of the respective trace 
operators, i.e. 

TR(G) :=7G[^i(|G|+i)], 
TR{D) ■.= jD[Hii\D\+i)]. 

Clearly the kernels of jg a-^d 7D are given by Hi{\G\ + i) and Hi{\D\ + i) respectively. This 
leads to the following definition. 

Definition 4.3 (Boundary data spaces). We define the boundary data spaces BD(G) and 
BD(L>) as 

BD(G) := FidGI +i)^«i(IGI+i) 

and 

BD(L') := Hi{\D\ + i)-Lffi(lo|+i). 
Lemma 4.4. The boundary data spaces are given by BD(G) = [{0}](1 — DG) and BD(D) = 

m]{i-GD). 

Proof. See |12l Lemma 5.1]. □ 
Theorem 4.5. The operators 

7g|bd(G) : BD(G) ^ TR(G) 

and 

7d|bd(D) :BD(I))^TR(D) 

are unitary. 

Proof. Let u G BD(G). Then we get for each v £ Hi{\D\ + i) that 

\{-fGu\v)\ = \ {Gu\v) + {u\Dv)\ 

= \ {Gu\v) + {DGu\Dv)\ 

= \{Gu\v)H,(\D\+i)\ 

< \Gu\Hi{\D\+i)\v\Hi{\D\+i) 
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4 Boundary Control Systems 



and hence, 

\lGu\H^l{\D\+i) < \Gu\Hi(^\D\+i) = \/|Gu|2 + \DGu\'^ = v^|Gu|2 + |m|2 = \u\H^(_\G\+i)- 

On the other hand we have 

{jGu\Gu) = {Gu\Gu) + {u\DGu) = {Gu\Gu) + {u\u) = \u\]j^^^Q\^^y 

which gives — l7G^^lH_i(|D|+i)- That 7g|bd(G) onto, follows by the definition of 

BD(G) and TR(G). The assertion for 7_d|bd(d) follows by interchanging the roles of D and 
G. □ 

Since G[BD(G)] C BD(L») and L'[BD(L>)] C BD(G) we may consider the following restrictions 
of G and D 



D : BD {D) BD (G) 

<p^ D(t> 

and 

G : BD (G) BD (D) 

(t>^G^. 

The operators D and G enjoy the following surprising property. 

Theorem 4.6. We have that ^ 

G] =D= Ig' 



In particular, G and D are unitary. 

Proof. See [12 Theorem 5.2.]. □ 

Remark 4.7. The operator 7_d|bd(D) G (7g|bd(G)) • TR(G) — )■ TK{D) is unitary and can 
be interpreted as an abstract version of the Dirichlet-to-Neumann operator. 

After these preparations we show, how systems with boundary control and boundary observa- 
tion can be treated within the framework of Section 3. For doing so, let C € L{Hi{\G\ +i); V) 
for some Hilbert space V and assume that F is given by 



F:=i^ Hi{\G\ + i) C Ho{\G\ + i) ^ Ho{\D\ + i) © K (4) 

As in the definition of an abstract linear control system Cm,f,b the adjoint of F comes into 
play. We compute it explicitly in the next theorem for the case, when G is assumed to be 
boundedly invertible. In applications this requirement can be guaranteed by assuming certain 
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geometric properties of the underlying domain (e.g. segment property, Lipschitz boundary 
and so on). 

Theorem 4.8. Let F be given as above and let G be boundedly invertible. Then 

F* : D{F*) C Ho{\D\ + i) © F ^ Ho{\G\ + i) 

where is the dual operator ofC with respect to the Gelfand-triplet Hi{\G\+\) C Hq[\G\+\) C 
H_i{\G\ +i) and 

D{F*) = {(C, w) G Ho{\D\ +i)eV\DC + C''w£ Ho{\G\ + i)}. 



Proof. See O Theorem 5.4.]. □ 

Remark 4.9. Let Cm,f,b be an abstract hnear control system, where F is given by (j4]). We 
assume that G is boundedly invertible. Note that, as a consequence, D is boundedly invert- 
ible as well. An element {x,((^,w)) G Hq{\G\) © (^Ho{\D\) (B belongs to the domain of 
-F* 



, if and only if x £ Hi{\G\) and + C^w € Hq{\D\). The latter is equivalent to 
r U J 

^DiC + iDr'C'w) = 0. (5) 
Recall from ([3]) that M is of the form 



M{z) 



for suitable operators Mi^y and K : B(c{r,r) L{Hq{\G\) © Hq{\D\) © V). Note that due 
to the block structure of F* , this operator has indeed four lines and columns. With (x,^) = 
(x, {C,,w)), the third and fourth line of the equation given by Cm.f,b read as 

do7rvK{dQ'^) (7r|^„(|G|)3; + '^Hn(im^) + doTTvK{dQ^)TryW + 7ryMi_i2y + Gx = vry/ + TryBu 




and 

Mi^2ox + Mi,2i7r^g(|^|)C + Mi^2iTTyW + Mi,22y = Try/ + vrySu, 
respectively. We may rewrite this as 

do7rvK{do^)TT*y 7ryMi,i2\ (w' 
Mi,2i7r^ Mi,22 J \ y 

Try/ + ttyBu - Mi^2ox - Mi,2i7r^^^|^|^C 



(6) 
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5 A Boundary Control Problem in Visco-elasticity 

or equivalently as 

( doTrvK{do^)7r*y -7rvB\ ( w\ 
V Mi,2i7r^ -ttyB J \ uJ 

_ / Try/ - 7rvMi,i2y - 9ovry/s:(9o"^)(7rHo(|G'|)^ + ^^o(|D|)''^ ~ 

If the material law M satisfies the solvability condition ([1]), then the operator on the left 
hand side of ([6]) is boundedly invertible and thus, we can express w in terms of x, / and u. 
Plugging this representation into ([5]) we obtain the boundary control equation. Analogously, 
by assuming that the operator on the left hand side of Equation ([7]) is boundedly invertible, 
we can express w in terms of x, f and y and hence. Equation ([5]) yields the boundary 
observation equation. 




5 A Boundary Control Problem in Visco-elasticity 

In this section we apply our results to a boundary control problem for the equations of visco- 
elasticity. For this purpose we introduce the required differential operators. Throughout let 
C M*^ be an open subset of M", n G N, n > 1. 

Definition 5.1. We denote by L2(f^)"'^" the Hilbert space of n x n-matrices with entries in 
L2(J^) endowed with the inner product 

($|^')i,(o)nxn := J trace (^>(x)*^(x)) dx G ^2(5^)"^") • 

Moreover let ffsyjn(J^) ^ ^2(^1)"^" denote the closed subspace of symmetric n x n matrices. 
We define the operator Grad as the closure of 



Grad|c-(n)n : C^{nr Q ^2(J^)" ^ ^sym(J^) 

'1 

and Div as the closure of 



i)ie{i,...,n} ^ (7: idj(t>i + di(t>j) 



i,je{l,...,n} 




where we denote by C^yjn(r2)"^" the space of symmetric n x n-matrices with entries in (7^(0). 
Furthermore we extend the meaning of grad by defining it as the closure of 

grad|c,c«(f,)„ : C,°^(0)" C ^2(0)" ^ L^inr""^ 
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and similarly div as the closure of 

lc-(fi) 



div|c-(mnxn : C,^(0)"^" C L2(f^r''" ^ L2{VlY 



«e{i,...,n} 



leaving it to the context to determine if the scalar or the matrix version of these operations 
are meant. 



An easy computation shows that Grad and Div are formally skew-adjoint, likewise the ex- 
tended operations grad and div are formally skew-adjoint. Following the notation intro- 
duced in Section 4 we define Grad := —(Div)*, Div := —(Grad)*, grad := —(div)* and 
div := —(grad)*. The equations of visco-elasticity are given by 

dlQx{t)-J}\wT{t) = f{t), (8) 

t 

T{t) = M Grad x(t) - j g{t - s) Grad x{s) ds. (9) 

— oo 



Here x G H^fi{M.\L2{^)'^) and T G H^^Q{K]Hsyra{^)) are the unknowns, denoting the dis- 
placement field and the stress tensor, respectively. The density function g G Loo(f^) is 
assumed to be real-valued and uniformly strictly positive, i.e. > ci > 0. The tensor 
M G L (i?syjn(r2)), linking the stress and the strain tensor, is assumed to be selfadjoint and 
satisfies M > C2 > 0. The function g : M>o — )■ C is assumed to be absolutely continuou^H, i.e. 
g{t) = Jq h{s) ds + go for some h G Li(M>o) and go G C. An easy computation shows that 
the convolution operator g* : H^^q(R; Hgyj^^Q)) — t- H,^^q(M; Hsymi^)) is continuous for each 

> with \\g * ||L(//,_o(M;H,ym(Q)) < h (|^Ili(R>o) + IS'ol) • Thus, for i/ > large enough, the 
operator (1 — M~^g*) is invertible, and hence we may write ([9]) as 

(M - g*)~^T = (1 - ]Vr^g*y^A'r^T = Gradx. (10) 



The boundary control and observation equations are given by 

TN = dox + V2u, 

TN = V2y-dox, (11) 
on dQ, where we denote by N the outer unit normal vector field. 

Remark 5.2. Since we have to compare Neumann- type traces and Dirichlet-type traces we have 
to determine a suitable control and observation space. For doing so, let us assume for the 
moment that dQ is smooth. We consider the space L2{dQ)^. We assume that the outer unit 
normal vector field N can be extended to 0, such that G Loo(f^;IK)" and divA^ G Loo(f^). 

•^The equations can also be studied in a more general setting, for instance g can attain values in L(_ffsym(f2)) 
(cf. [15]). 
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5 A Boundary Control Problem in Visco-elasticity 

For f,g & BD(grad) we formally compute using the divergence theoreixQ 
J f-gdS = J {f-g){N.N) dS 

3n dQ 

= 11 k,9*) -NdS + ^l [{Nkg*)k, f)-NdS 

an dQ 

= i 1 div mm k,9) + ^ / {{Nk9^)k, r) 

div (/iA^fe) . ^)l2(Q)" + ^(•^»^fc)^fcls^^^f)L2(n)"X" 



1 



2 1 \^^BD(div) {fi^k)i^k 



grad g) + ( grad / 

BD(div) 



BD(div) , 



This leads to the following choice for the control space U. 

Definition 5.3. Let N G Loo(J^;IK)" be such that divA^ G L^o- We define the bounded linear 
operator u : i?D(grad) — t- BD(div) by uf := vrBj^^div) (/i-^fc)j fcgji „}■ We assume that the 

operator div i' + i'* grad is positive, i.e. for every / G BD(grad) \ {0} we have 

• \ 

div V + V* grad I / 



BD(grad) 



We define the Hilbert space U as the completion of BD (grad) with respect to the inner product 
(•!•){/ : BD(grad) x BD(grad) ^ C 

/ BD(div) / 



{f,9)^l{{iyf 



grad g ) + ( grad / 

BD(div) 



We denote the embedding BD(grad) [/ by t. 

In the following we require that Korn's inequality holds, i.e., Hi{\ Grad | +i) ^ -f^id grad | +i) 
(for sufficient criteria see [1] and the references therein). We consider the bounded operator 
j : BD(Grad) U given by j = t o 7rBD(grad) ^BD(Grad) ^'^d compute 

{jf\9)u = (^BD(grad)«^7rBj3((.j.^^)/|g)c/ 



2 I ( '^^BD{grad)«7rBj3(Grad)/ 



grad g) + { grad ^BD(grad)K^BD(Grad)/ 



1^9 



BD(div)^ 



/ 



7rBD(Grad)'^*7rBD(grad)^* grad fil 



BD(Grad) 



*Note that due to the assumptions on the vector field A'', the matrix-valued function {fiNf;)i^ks{i,...,n} hes in 
_D(div) for each / G Z)(grad). 
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+ / 



7rBD(Grad)/«*7rBD(grad) ^iv l^g 



I BD(Grad)^ 



/ 



^VrBD(Grad)«^*Vr*j3(gi.ad) ( ^* gl'ad + div 1 fif 



BD(Grad) 

for each / G BD(Grad), g € BD(grad). This gives 



^^BD(Grad)'^*^BD(grad) g^ad + div 



and, consequently, 



7DivVrBD(Div) Grad fg (12) 
= ^7DivVrBD{Div) Grad vrBD(Grad)K*7r^D(grad) grad + div i/^ g. 

Remark 5.4. We give a possible interpretation of the latter equality. For this purpose we 
compute formally using the divergence theorem 

(Gradr5)iV-/d5 



do. 



^ (7DivVrBD(Div) Grad j*5f|7r^D(G^^d)/) 

■ 2(7DivVrBD(Div) Grad VrBD(Grad)/«*7rBD(grad)'^* grad C/|7rBD(Grad)/) 
+ 2(TDiv7rBD(Div) Grad 7rBD(Grad)/«*7rBD(grad) ^iv J^^kBDCGrad) /) 

1 • 1 • • • 

^ 2^^BD(Grad)«^*7rBD(grad)^* g^ad g|/) + -(Grad vtbd (Grad) '^*7rBD (grad) g^ad g\ Grad /) 

1 • 1 • • • 

+ 2^^BD(Grad)«^*vr*j3(grad) ^iv vg\j) + -(Grad vrBD(Grad)«*7r*D(g^^j) div vg\ Grad /) 

■ 2^^BD(Grad)«^*^BD(grad)^* gl'ad 9|/)BD(Grad) + 2 ^^BD(Grad) '^*7rBD(grad) ^iv I^fi'l/)BD(Grad) 

1 • ^ 1 • 

: -(grad 5l'^7rBD(grad)'^7rBD(Grad)/)BD(div) + 2(^5! grad 7rBD(grad)'^7rBD(Gj.ad)/)BD(div) 

■ (5'KBD(grad)'^7rBD(Grad)/)c/ 



an 



9 ■ 7rBD(grad)'^7rBD(Grad)/ ^5" 
J 9 ■ 7I"BD(grad)'^7rBD(Grad)/ ^5" + j 9 ■ {I - T^BD{gri,d))t^T^BTi{Gvad)f 

an dvi 

= y fir - KTTbd (Grad)/ ^5" 

an 

= j 9-fdS 
an 
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for each / E BD(Grad), (7 G BD(grad). Hence, equality ()12p can be seen as a generalization of 

(Grad j* g) N = g on dQ (13) 
to the case of non-smooth boundaries. 

We now want to transform the equations ([8]), (|10p and ()lip into a system of the form treated 
in the previous subsections. In the terminology of Subsection 4 the operator Grad should play 
the role of G and Div the role of D. Since we have assumed in Theorem 14. 81 that G is boundedly 
invert ible, we require that Grad[L2(0)"] is closed in i?sym(f^)ll 

Then the projection theorem 

yields the following orthogonal decomposition^ 



L2inr = [{0}] Grad ©Div [H,y^{n)], 
H,yn,{n) = [{0}] Div ©Grad[L2(l^)"]. 



We define the orthogonal projections VTDiv : -^2(^)" Div [-ffsym(^)] and vrQi-ad ■ -ffsym(f^) — ^ 
Grad[L2(r2)"']. Note that due to the closed graph theorem the operator Grad := vrcrad Grad vr^j^ 
is boundedly invertible and so is (Grad) = — vroiv Div t^q^.^^- Furthermore let us denote by 



'-Grad the embedding Hi 
ary problem 



Grad 



do 



( VTDiv £'7rJ,i^ ( 

0\ /vrGrad(M-5*)- 

V (0 



+ i ) ^ ¥l\ (|Grad| + i) . We consider the following evolution- 



\ ( 

+ 



V 







^/2) 



+ 





^ ) 





Gradj 




0) 



-C" 












\ 



\ 



/ / \ 











-1 



(14) 



where C : H\ 
given by 



Grad -|- i j — )■ ?7 is given by Cx := i'7rBD(Qi.ad)''Grad3^- The material law K 

Kiz) = ( TTGrad (M - V2^g{-iz~')y^ vr^^^^ 
z 



IS 



^The closedness of the range Grad[L2(n)"] holds, for instance if Hi{\ Grad | + i) is compactly embedded in 
1/2 (i^)" and we refer to for sufficient conditions on yielding this compact embedding. Note that this 
compact embedding yields then a Poincare-type estimate, which in turn yields the closedness of the range. 
Thus, the minimal assumption is the validity of a Poincare-type estimate. 

o 

®Note that the closedness of the range of Grad also yields the closedness of range of Div. Since we do not 
want to give the details of the proof here, we use the closure bar for convenience. 



18 



and satisfies the solvability condition ([T]). Indeed, using the representation z ^ = it + for 
some V > ^jteMifzG Bc{r, r) we estimate 



and 



= IHez-Vcrad ( E (2^)^ M-^5(-iz-i)M M-Va,ad 

Vfc=0 / 
= ^^TTGradM-^Grad 

+ 9^e^Grad^-' (^/2^M-i?(-i^-^)) (f;(27r)Uf-^^5(-iz-i)M M-V^..^ 



= Z^TTGradM V^j-ad 

ll^-'f (|/i|Li(R>o) + Iffol) 
>Z/C2 ^ " ^ 



Vfc=0 



l-l/-l||Af-l|| (|/l|L,(R>o) + lffO 



where we have used g{ — 'iz ^) = zh{—\z ^) + -^gQ. Summarizing this gives yit z ^K{z) > 1. 

The operator J is given by J = / \ and thus, || J|| = Since iHeMi^22 = 1 Theorem 
13.21 applies and thus, the control system given by (|14p is well-posed. Next, we compute C^. 



For that purpose let x £ H 



Grad 



+ i ) and w ^U. Then 



{C''w\x) = {w\Cx)u 

= («^li7rBD(Grad)''Grada;)(7 

= (^BD (Grad)i* ^ I ''Grad a:) ( | Grad | +i) 

= (^BD(Grad)J*"'kGradx) + {GxiidTT*^j^(^^^^^f w\ Grad iGradx) 



Gradx j 



(^BD(Grad)J*^|2;) + (^7rGradVrBD(Div) Grad J* U- 

(^BD(Grad)J*^|a;) + ( (Grad) ^GradvrBD(Div) Grad J* u; 



Summarizing, we get that = 7rBD{Grad)-^* ^ ^Grad) '''"Grad7rBD(£)iv) Grad j*. According to 
the definition of the domain oi ( — (^Grad^ —C^\ the implicit boundary condition for the 
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5 A Boundary Control Problem in Visco-elasticity 



system reads as 



(G^dY T + Cw e Ho 



Grad 



+ i . 



Hence, 



(Cvad) 



+ i]OHi 



Div 



+ i C i^i (|Div| + i) 



From 



7rGradvrBD(Div) Grad j*w = vrcrad Grad7rBO((jj.^^)j*u; G i7i(| Div | + i) 
it thus follows that T £ Hi (| Div | + i) and 

-1 



fCrad) ) C'w 
1 



iDivT = 7Di 

= 7Div ( - ((^Gradj j 7rGradvrBD{Grad)i*^ - ^GradvrBD(Div) Grad f 



w 



= -7Div7rGrad7rBD(Div) Grad j*'u; 
= -7DivVrBD(Div) Grad fw, 

o 

where we used that 7Div vanishes on the domain of Div, which is a superset of the domain of 



(^Gradj . Using ([6]) we get 



1 

V2 1 



w 



-V2u-Cv 



and hence 



w = u 



Cv = -V2t 



Analogously one obtains, using ([7]), that 



1 V2\ fw 
I )\u 



J^BD(Grad)''Grad'y- 



-Cv 

-y 



and thus. 



w 



Cv - V2y = j7rBD(Gr 



id) '•Grad 



V2y. 



Following the reasoning of Remark I5.4| we may interpret the resulting boundary control and 
observation equations as 

TN = v-u = dox + V2u, 
TN = V2y -V = \f2y - Oqx 



on 50. 
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